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Abstract 

Reactions at a strongly segregated interface of a symmetric binary 
polymer blend are investigated via Monte Carlo simulations. End func- 
tionalized homopolymers of different species interact at the interface in- 
stantaneously and irreversibly to form diblock copolymers. The simula- 
tions, in the framework of the bond fluctuation model, determine the time 
dependence of the copolymer production in the initial and intermediate 
time regime for small reactant concentration poR^ ~ 0.163 • • ■ 0.0406. The 
results are compared to recent theories and simulation data of a simple 
reaction diffusion model. For the reactant concentration accessible in the 
simulation, no linear growth of the copolymer density is found in the initial 
regime, and a \/f-law is observed in the intermediate stage. 



1 Introduction 

Reactions at interfaces have attracted great interest because of a variety of prac- 
tical applications. They are involved in facilitated transport across biological 
membranes, and in selective removal of a component of a mixture by contact 
with a reactive liquid Recent interest stems from extensive commercial ap- 
plication of reactive polymer blending technique || and progress of the ex- 
perimental studies j|, |, 1^. Melt blending of different polymer species (denoted 
A and B) is widely employedQ to improve the application properties of poly- 
meric materials. Since the entropy of mixing of high-molecular-weight polymers 
is extremely low, macromolecular alloys can often be described as an assembly 
of interfaces between relatively pure phases. Block copolymers containing both 
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types of monomers can be used to tailor the properties of polymer-polymer in- 
terfaces. They promote mixing by means of reducing the interfacial tension and 
the droplet-droplet coalescence rateQ and mechanically strengthening the com- 
posite material due to entanglements across the interface. An effective method 
of creating copolymers at the interface is in-situ production[||, e.g. chemical 
reactions of end functionalizcd homopolymers. 

Recently, reactions at a flat interface in a binary, symmetric polymer blend 
have been studied independently and in some detail by O'Shaughnessy and 
Sawhneyl^, ^ and Fredrickson[0 and Milner[|llj: A and B polymers are as- 
sumed to be immiscible, but structurally symmetric. The almost pure phases 
A and B host a small fraction po of chains with a terminal reactive group. 
The A reactive ends can only react with B reactive ends and form AB diblock- 
copolymers at the interface between the bulk phases. The reactions are instan- 
taneous and irreversible. The interfacial density a of copolymers at the interface 
grows according to the rate equation: 

% - m,i (1) 

where t measures the duration of the reaction and po denotes the time inde- 
pendent fraction of reactive chains in the bulk. The reaction rate K is largely 
independent from the details of the microscopic mechanism, because subdiffu- 
sive monomer motion on small time scales guarantees reaction whenever the 
coils of two reactive chains over lap ||T^. For small number density po of reac- 
tive chains {poRg 1 with Rg radius of gyration) and high incompatibility, 
Fredrickson and Milner|jll| predict three time regimes: At times smaller than 
T = D/KqPq, the density of reactive chains at the interface remains close to its 
bulk value. The reaction rate within the framework of the Rouse model is 
time- independent, and given by||, |lo): Kq « 100.61? i?^/ In iV. It follows that 
the interfacial density of copolymers grows linearly with time in this regime. At 
times larger than t[|i^, there is a depletion hole of reactive chains with spatial 
extension V Dt, and the reaction rate is determined by the flux of reactive chains 
to the interface. The reaction rate decays as K ~ y^D/irpQi, i.e. the interface 
acts much like an absorbing boundary. The copolymer density grows like t^^"^ 
in this regime. For even longer times t > ^^Rg/N^Dpg ($ monomer density), 
the copolymer accumulates at the interface and forms a brush. This prevents 
reactive chains from reaching the interface, and the reaction rate becomes ex- 
tremely small, decaying as 1/t^/lnt. The copolymer density increases only as 
Villi. 

In spite of the practical applications of interfacial reactions, qualitative ex- 
periments to determine the time dependence of the copolymer concentration 
at the interface are still rare|^, ||, ^, |l5j. Their interpretation is complicated 
by structural asymmetries between the two homopolymer species, or the reac- 
tive and inert polymers of the same species. Furthermore, large differences in 
the mobilities of the species lead to qualitative modifications of the reaction 
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kinetics [p^. In such a situation, Monte Carlo simulations can provide valuable 
information for understanding the reactions in these inhomogeneous, complex 
fluids and for exploring the validity of approximations in the analytic treatment. 
Such simulations of reactions at surfaces and of catalysis of small molecules have 
proven extremely useful in studying the stochastic dynamics of reaction-diffusion 



processes 17 . The present simulations are the first to investigate reactions at 
polymer interfaces, and they cover the initial and intermediate time regime. 
They yield the time evolution, as well as the spatial depletion profile, of reac- 
tants at the interface. We compare our results to the theories described above 
and to a simple reaction-diffusion model. 



2 Model and simulation technique 

If the microscopic reaction is sufficiently fast, the reaction rate is largely in- 
dependent of the microscopic mechanisms [|[ |l^. Since the study of reaction 
kinetics poses high computational demands, we employ a coarse grained and 
highly computational efficient lattice model of a dense polymeric melt. The 
Bond Fluctuation Model (BFM)||l^] retains the relevant features of polymeric 
materials: connectivity of the monomers along a chain molecule, excluded vol- 
ume interaction of the monomeric units, and a short range thermal interac- 
tion potential. The BFM has proven very useful for studying various related 
properties of dense polymeric melts. Among them are single chain dynam- 
ics in athermal melts p9|, pO| and binary blends [pT|, living polymerization , 
adsorption in polymer brushes ||2^, and interfacial properties in binary ||2^, 
and ternary blends |^^. In the framework of the BFM each segment occupies 
all eight sites at the corners of a unit cell in a simple cubic lattice, and no 
site can be doubly occupied. Adjacent monomers along a chain molecule are 
connected via one of 108 bond vectors. These are chosen such that the local self- 
avoidance prevents chains from crossing each other during their motion. The 
fluctuating bond length permits the implementation of random, local monomer 
displacements [^Sj such that the dynamics of an isolated chain is well describable 
by the Rouse model[^ We work at a filling fraction 8$ = 0.5 of occupied 
lattice sites, and all length are measured in units of the lattice spacing u. At 
this density the BFM reproduces many characteristics of a dense polymer melt. 
The blend consists of two species of homopolymers, denoted A and B, of equal 
chain length N = 32, and diblock copolymers AB of chain length 64. This is 
well below the chain length N = 200, where first signs of repation-like motion 
are observed. We employ a square well potential, which encompasses the first 
peak of the pair correlation function, i.e. comprising all 54 neighbor sites up to 
a distance -\/6 lattice units: e — cab = —^aa = —£bb = O-lksT. Monomers of 
the same kind attract each other, whereas there is a repulsion between unlike 
species. This incompatibility corresponds to ~ 17, and is well inside the 
strong segregation limit for binary blends. Both, the interfacial properties of 
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the binary blend ||2^, ^ as well as the equilibrium behavior of copolymers at 
interfaces ^ have been investigated extensively. 

We work in a 64 x 64 x 128 geometry, and for the highest fraction of reactive 
chains additional data were obtained for system size 64 x 64 x 256. Periodic 
boundary conditions in all three spatial directions were applied, such that the 
systems contain two interfaces parallel to the xy-plane. The initial configura- 
tions for the binary blend have been carefully equilibrated using 10^ local MCS. 
(On the average each monomer attempts to jump once within a Monte Carlo 
Step/MCS). After this initial equilibration of the interfacial properties, the equi- 
libration was continued and starting configurations for the reaction study were 
collected every 50,000 MCS. For the smaller system size we average our results 
over 96 independent simulation runs, corresponding to a total of more than 3 
million monomers. For the larger size we employ 30 independent configura- 
tions. Simulations were performed on a CONVEX SPP1200, using a trivial 
parallelization strategy. 

Some interfacial properties of the starting configurations are presented in 
Fig. I. The blend is well segregated, i.e. the coexisting phases consist almost 
entirely of A or B polymers. Between the coexisting phases there is an inter- 
face of width w ~ 3.7. The conformational properties have been discussed in 



detail in Ref.|2J, 25 . We simply emphasize that, in agreement with mean field 
theories 25 , there is a pronounced segregation of chain ends at the interface 
and a concomitant depletion at a distance Rg. The total monomer density is 
somewhat reduced at the center of the interface. 

In these equilibrated starting configurations a fraction of homopolymers was 
randomly chosen to bear a reactive end, and the reaction was started at time 
t = 0. This corresponds to an experimental situation in which a (fast diffusive) 
reactive agent/catalyst[^ is added to the blend, or the reaction is initiated by 
radiation. O'Shaughnessy and Sawhney^, ^ analyzed the qualitative depen- 
dence on the monomeric reaction probability. For very low monomeric reaction 
probability Q, correlation in the density of reactive particles are only weakly 
perturbed by the reactions, and the reaction rate K is proportional to the equi- 
librium contact probability of reactive ends. However, this mean field kinetics 
is only appropriate if the reaction probability is smaller than a crossover value 
Q* ^ l/lnA^m ||. Therefore, reactions in high molecular weight polymers 
are typically diffusion controlled. This diffusion controlled reaction kinetics is 
largely independent from the details of the monomeric reaction mechanism, es- 
pecially from the monomeric capture radius and reaction probabilities. Since 
the simulation can treat only rather short polymers, we use a large monomeric 
reaction rate Q and capture radius to make sure that we work in the controlled 
regime. If the distance between reactive A-ends and B-ends is less than the 
capture radius ^6, reactive homopolymers form AB diblock copolymers instan- 
taneously and irreversibly. Of course, the concrete value of the capture radius 
•\/6 has no physical significance; it is the largest value the model allows to be 
compatible with local monomer dynamics and chain connectivity. According to 
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analytical treatments l^], ^ |T^, |ri|, we expect our results to be characteristic for 
diffusion controlled kinetics and rather independent of the concrete monomeric 
reaction mechanism (e.g. the model also applies to long polymers with a large 
but finite monomeric reaction rate). 

Otherwise the reactive chains behave the same as the non-reactive polymers. 
We considered systems with three different values of the initial number density 
of reactive chains in their corresponding bulk phases: po — 1/2048, 1/4096, and 
1/8192. With the measured radius of gyration Rg = 6.93, the scaled number 
density ranged in our simulation from poRg = 0.041 to 0.1625. Furthermore, 
we determined the self-diffusion constant from the mean square displacements 
in the plane parallel to the interface for the largest system size to = 1.17(5) • 
10~^u^/MCS, which is very close to the corresponding athermal value. 



3 Simulation results 

Upon activating the reactive polymers, AB diblock copolymers are created at the 
interface. The time dependence of the interfacial density for different fraction 
of reactive chains is presented in Fig. ^ For the highest fraction, the simulation 
data for the two system sizes agree to within the statistical accuracy. Since 
finite size effects are expected to be even smaller for lower fractions of reactive 
chains, they can be neglected for our data. In contrast to the predictions above, 
the simulation data for the copolymer interfacial density do not exhibit a linear 
growth law at small times. Due to the high initial reaction rate in the Monte 
Carlo simulation, the theory underestimates the interfacial density over the 
entire time covered by the simulation for the two largest reactant concentrations. 
The broken lines in the figure show the theoretical prediction for the diffusive 
growth in the intermediate time regime. For the lowest density, the simulation 
data are rather well described by the prediction for the intermediate regime. 
However the agreement is fortuitous, as the simulation data cover mainly the 
initial time regime (t w 3 • 10^ MCS). 

To eliminate the influence of the fast initial rate, it is instructive to monitor 
the growth rate K{t). The solid lines in Fig. ^ show fits to the time development 
of the copolymer density which are used to determine the reaction rate. These 
rates are presented in Fig. ^ in scaled form. The solid lines show the simulational 
results, while the dashed ones correspond to the theoretical predictions in the 
initial {t < r) and intermediate {t > r) regimes. First order corrections 11 to 



describe the crossover between the two regimes are also included into the figure: 

— — ^ = 1 1= — tor t <C T mitial regime (2) 

da/dt 1 ftY^'^ ftV""^ r , . ^ 

— ' ' 0.254 — 1 for t ^ T intermediate regipi^ 
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Even for poR^ as small as 0.0406 the prediction underestimates the reaction rate 
by an order of magnitude, and there is no sign of a constant reaction rate at small 
times. Furthermore the simulation data for different fractions of reactive chains 
do not collapse onto a single curve. In the intermediate time regime, however, 
the simulation data are consistent with the theoretical prediction K ^ 
and the anticipated proportionality constant is in rough agreement with the 
simulation data, though the asymptotic limit t/r 3> 1 is not reached for the 
lower fractions of reactive chains. 

In Fig. ^ we present the non-equilibrium profiles across the interface for 
PoRl — 0.163 and t/r fa 24.4. The copolymers accumulate at the interface, 
and their composition profile resembles its equilibrium shape, as determined in 
Ref. The copolymers do not yet form a brush at the interface. Therefore 
the simulation data are still in the intermediate regime. The homopolymer 
density profile, however, is altered as they are displaced from the interface. The 
simulation data also confirm the predicted depletion of reactive chains in the 
interfacial zone, and the growth of this depletion hole with time. Thus the 
simulation data corroborate the main features of the theory in the intermediate 
time regime, but disagree strongly with predictions in the initial stage. 

There are several possible reasons for the difference in the initial stage be- 
tween the simulation data and the theory: 

The local equilibrium structure of the homopolymer interface exhibits an 
enrichment of chain ends at the interface. Therefore the initial reaction rate 
is expected to be higher than that for a uniform initial distribution of reactive 
segments. This effect has been used to rationalize the enhanced reaction rate in 
PS-PMMA interfaces |T^. However, it is not clear to what extent this correlation 
in the initial condition accounts quantitatively for the discrepancies between 
simulation and theory. Futhermore, the theory deals with Rouse chains at an 
interface in the Helfand regime 1 ^ x ^ This latter condition can be 

satisfied only for very long chains, whereas the chain length N must not be too 
large so that the chain dynamics is describable by the Rouse model. 

Additionally, the reactant density PoRg — 0.0406 in the simulations might 
still not be small enough to meet the condition PoR^ <C 1 assumed in the 
theory. The constant reaction rate should be observable for times much larger 
than the relaxation time of the internal chain motion, but much smaller 
than the crossover time t to the intermediate regime. Note however, that the 
crossover time scale r and the Rouse time = GRg/n^D w T{78.QpoRg/ \nN)^ 
are not well separated for the reactant concentration studied. Unfortunately, 
simulation at even lower fraction of reactive chains pose very high statistical 
demands, which are beyond our computational facilities. 
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4 Reaction diffusion model 



We decided to compare our results to additional simulations of a simple random 
walk model. Bulk reactions [p^ in polymeric systems can be modeled schemat- 
ically by random walkers which react if their mutual distance is smaller than 
the radius of gyration in the corresponding polymer system. Therefore we con- 
sider two types of random walkers, denoted A and B, on a simple cubic lattice. 
A-type walkers are restricted to the half-space z > 0, whereas B-type walkers 
are confined to z < 0. An A-type walker annihilates with a B-type walker 
(and vice versa), whenever their mutual distance is smaller than R. The ini- 
tial time regime for a class of similar models has been discussed by Burning 
and O'Shaugnessy fl^ . The model is an appropriate caricature of polymer 
reactions at interfaces for the initial and intermediate time regime. Since it 
eliminates the shortest time scales, which are associated with the segmental 
motion on distances smaller than the radius of gyration, and also neglects the 
inert homopolymer matrix, much lower volume fractions for reactive particles 
are accessible in the simulations. Furthermore, many of the above-mentioned 
complications characteristic of polymeric systems are absent. We work with a 
capture radius i? = 8 and a 256 x 256 x 512 geometry, and a 1024 x 1024 x 2048 
system for the lowest density. Initially, the random walkers are uniformly dis- 
tributed in their volumes. We employ at least 512 independent systems. 

The reaction rate in its scaled form as well as the prediction of Fredrick- 
son and Milnerjl^ are presented in Fig. |5[ For very small density of walkers 
PqR^ < 0.002, the simulations show an almost constant initial reaction rate, 
in agreement with the predictions. The crossover to the intermediate regime, 
and the diffusive growth in the intermediate regime, is described by the theory 
in the limit of vanishing reactant concentration, and the corresponding scaling 
function is obtained as the lower envelope of the simulation data. The spatial 
density profile (cf. Fig. |5| inset) in the intermediate time regime exhibits a pro- 
nounced depletion hole, and is qualitatively similar to the profile obtained in 
the simulation of the polymer system. In the vicinity of the interface there are 
some deviations due to the finite interfacial width in the polymer system. 

For concentrations comparable to those used in the simulation of the poly- 
mer systems, the reaction diffusion model shows similar deviations, i.e. a de- 
caying reaction rate in the initial stage which is always larger than in the limit 
Pq — > 0. However the deviations are somewhat smaller than for the polymer sys- 
tem. This is partially due to the fact, that the reaction in the polymer system 
does not take place instantaneous when the chain volumes of the reactive pair 
overlap, but the reaction is retarded by a time of the order r^. Therefore we 
conclude that the finite concentration of reactive chains accounts qualitatively 
for the discrepancies between the analytic theory and simulations of the poly- 
meric system. However, the simple reaction-diffusion model cannot describe the 
polymeric system quantitatively for finite reactant concentrations. 
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5 Conclusions and outlook 



In summary, we have presented extensive simulations of reactions at interfaces, 
both for polymeric systems as well as for a simple reaction-diffusion model. We 
have compared our simulation results to recent theoretical predictions [|[ |[ ^ 
0. At intermediate times wc find agreement between the theory and the sim- 
ulations of both models. The initial state and the crossover to the intermediate 
regime is described by the theory only for extremely small concentration of re- 
active particles. These are not accessible in the simulation of the polymeric 
systems, but are in the simulation of the random walker model. Results of the 
latter are compatible with the theoretical predictions in the limit poR^ — s- 0, 
and the simulations yield estimates for the complete crossover scaling behav- 
ior of the reaction rate. However, for non-zero but small concentrations of 
reactive particles, the reaction rate is grossly underestimated by the theory. 
Even for the smallest fraction of reactive chains simulated in the present study, 
PoRg = 0.0406 , the initial reaction rate obtained from the simulation is about 
an order of magnitude larger than predicted, and is not constant, but decays 
with time. It should be emphasized that our results are not an artifact of the 
small chain length used in the simulation. The reactant concentrations corre- 
spond to weight fractions of 25%, 12.5% and 6.25%. These weight fractions are 
within the experimentally relevant range. If one increases the chain length at 
fixed weight fraction, then poRg increases as \/N taking the initial condition 
even further from the assumptions of the theory. Therefore we expect that ex- 
periments will also observe a decay of the reaction rate over the entire time 
regime, unless extremely small volume fractions of reactive chains are used. 

Further studies might investigate the crossover between the diffusive growth 
and the late regime, in which the copolymer brush hinders the reactive chains 
to reach the surface. This regime is accessible for high concentration of reactive 
chains. At high incompatibilities the copolymers cannot reduce the surface 
tension completely, and there is a first order transition to a lamellar phase at 
high copolymer concentration [p6| . The direct measurement of the interactions 
between copolymer "coated" interfaces [|6| is an interesting topic, because it 
relates to the important effect of a reduced droplet-droplet coalescence rate. 
In the weak segregation limit, however, the addition of copolymer causes a 
continuous transition to the disordered phase to be encountered, (the Scott-line 
which exists for 2 < < 3 within mean field theory [^). Thus it would 
be interesting to study the temperature dependence of the reaction kinetics at 
weak and intermediate segregations, where reactions might lead to complete 
miscibility. 
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1.2 




Figure 1: Interfacial properties of the equilibrated starting configuration: The sohd lines present the 
normed density of A and B monomers, wheras the symbols denote the total monomer density (circles) 
and the density of chain ends (squares) scaled by a factor N/2. Note the enrichment of chain ends at the 
center of the interface. 
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0.008 




t in MCS 

Figure 2: Time dependence of the copolymer concentration at the interface: circles denote simulation 
data for po ~ 1/2048 and 64 x 64 x 256 system size, squares correspond to po = 1/2048, triangles to 
po = 1/4096, diamonds to po = 1/8192 and system size 64 x 64 x 128 respectively. The solid lines are 
fits to the simulation raw data, whereas the dashed lines represent the predicted time development in the 
intermediate regime. 
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Figure 3: Scaled reaction rate for the polymer system: The solid lines are the results of the 

Monte Carlo simulation for different values po-R^ as indicated in the figure. The two set of dashed lines 
correspond to the predictions of Predrickson and Milner[ll] for po-Rg 1- The higher values correspond 
to the asymptotic results for small and large values oft/r; the lower ones present first order corrections. 
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Figure 4: Non-equilibrium density profiles for poR^ — 0.163 and t/r ~ 24.4. The solid lines present the 
monomer density of the A and B segments of the homopolymers. The dashed lines correspond to the A 
and B monomers of the copolymers, whereas the monomer density of the two types of reactive chains 
is denoted by circles. The inset shows the ratio of the monomer density of reactive chains to the total 
monomer density of homopolymers. 
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Figure 5: Scaled reaction rate for the reaction diffusion system: The soHd hnes are the result of the Monte 

Carlo simulations, the values of poRg being indicated in the figure. The dashed lines arc the same as in 
Fig. 3. Only for very small concentration does one observe a reaction rate which is in agreement with the 
predictions [9]. The initial reaction rate is Kq w 157 

The inset presents the density profiles of reactants for poRg = 0.0625 and t/r = 0.28,1.13,4.53,18.12 
and 72.49. 
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